1 



Time of arrival in Classical and Quantum Mechanics 



Juan Leon 



Instituto de Matematicas y Fisica Fundamental, CSIC 
Serrano 1 13-bis, E-28006 MADRID, Spain 
leon @ imaff .cfmac .csic.es 



February 1,2008 



Abstract 

The time of arrival at an arbitrary position in configuration space can be given 
as a function of the phase space variables for the Liouville integrable systems of 
classical mechanics, but only for them. We review the Jacobi-Lie transformation 
that explicitly implements this function of positions and momenta. We then dis- 
cuss the recently developed quantum formalism for the time of arrival. We first 
analyze the case of free particles in one and three space dimensions. Then, we 
apply the quantum version of the Jacobi-Lie transformation to work out the time 
of arrival operator in the presence of interactions. We discuss the formalism and its 
interpretation. We finish by disclosing the presence (absence) of "instantaneous" 
tunneling for thin (thick) barriers. 

PACS: 03.65.Bz, 03.65.Ca, 03.65.Nk 

Keywords: time; phase space; Hilbert space; positive operator valued measures; 
tunneling. 



Classical and Quantum Mechanics use the notion of Newtonian time, an universal 
parameter that rules the evolution of all the dynamical systems of the universe. New- 
tonian time is "a priori" external to everything, physical systems and observers alike. 
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However, in many instances there are true time-like properties in the physical systems 
under study. In general, the answers to questions like: How long will it take to ...? 
or. When will it ...? etc. come in the form of a time that genuinely depends of the 
very system. The crux of the matter is that of finding the time (the time elapsed, or 
the instant in time) in which some property of the system will take a specific value, 
something that could be generically termed as "the time of arrival at that value". In the 
next section we deal with the formulation of this question in classical dynamics. The 
much more involved case of translating a time parameter into an operator on the Hilbert 
space, as requires the quantum treatment, is worked out in sections 2 and 3 for free and 
interacting particles respectively. In section 4 we point out some eccentric properties 
of the time of arrival at places that are classically forbidden. 



1 Deriving time in phase space 

The treatment of time as a phase space variable is a time-honored procedure. The term 
extended phase space was coined for the approach in which, to the n pairs {q,p) of 
the phase space variables of mechanical systems with n degrees of freedom, one adds 
the additional conjugate pair {t,pt), which requires the constraint pt + H{q,p) = 
for consistency. It seems possible to dismantle this construction trading in the way 
the pair {t,pt) by another canonical pair {q,p). Naively, one would single out one 
of the phase space variables {qi for instance) and make it equal to some parametric 
value (i.e. qi = x). Then, its canonical conjugate (the momentum pi in this example), 
would be fixed by the constraint, giving pi — (f){x\ q2,P2, ■ ■ ■ ,Pt)- The phase space 
would now be given in terms of {(92,^2), ■ • ■ , (9n,Pn), {t,Pt)}, with x acting as an 
external evolution-like parameter Hence, t{x) or in words, "the time of arrival at x" 
would be a legitimate question to ask. In spite of its apparent generality, is seldom 
possible to accomplish this program, not because of its very difficulty, but due to the 
non-fulfillment of some of the many conditions necessary for the existence of solutions. 
Here we will discuss the case of integrable systems for which there is a global construct 
for t{x), that we will describe explicitly. 

In the modern approach to classical dynamics (a standard reference is [jl|], a very 
readable text can be found in |Q]), a Hamiltonian system is called completely integrable 
(a la Liouville) when it satisfies the conditions a and b below; 

a. There are n compatible conservation laws 

<I>i((ji, . . . ,9„,pi . . . ,Pn;t) = Ci, « = 1, . . . ,n, that is: 

a.l. = {$„H} + ^ =0, V^ = l,...,7^. 

a. 2. = 0, Vz,j = l,...,n. 

b. The conservation laws define n isolating integrals that can be written as: 

b. l. $j = Ci ^ Pi = (f>iiqi, . . . Ci, . . . ,C„;t), V « = 1, . . . ,n. 
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In these conditions, the solution to the Hamihon equations is an integrable flow, de- 
scribed by a system of holonomic coordinates {q{t),p{t)) in phase space for each in- 
stant of time: 

qi{t) = qi{qo,Po;t), i = l,...,n. (1) 
Pi{t) =Pi{qo,Po;t), i = l,...,n. 

In particular, given the set of initial conditions (go ,Po), the system arrives at the point 
{q{t),p{t)) of the phase space in the path independent instant t. Conversely, these 
points define the corresponding times of arrival. In this case, time meets the require- 
ments to qualify as a derived variable in phase space, whose explicit construction oc- 
cupies the rest of this section. 

For integrable flows there is a special choice of phase space coordinates that math- 
ematically eliminates the effects of interactions (because the new positions are ignor- 
able coordinates). In other words, integrable systems are canonically equivalent to a 
set of translations (or of circular motions) at constant speed. It is customary to de- 
note the variables that determine these translations as action-angle variables, a name 
which strictly is appropriate only for the case of periodic systems, where the (closed) 
flow lines are topologically equivalent to circles. For these integrable flows, there is a 
canonical transformation W (the Jacobi-Lie transformation) to free-Uke variables 

{q,p; H{q,p)} ^ {Q, P; Ho{P)} (2) 

where H{q,p) = Hq{P). The most useful form of this transformation is P), that 
is, a function of the old positions and the new momenta, so that 

dW{q,P) dW{q,P) . 
Qi = — OP — , Pi = 1 = 1,. ..,n (3) 

The choice Hq{P) = 2_,j 2m relates the free coordinates Pi{t) = Pi and Qi(t) = 
+ Qi of the translation flow to the positions and momenta {qi{t) , Pi{t)) of the 
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actual flow generated by = J2i 1^ + ^(9)- In this work we shall only consider 

unbound systems with positive energy H = Hq > 0. For this reason we choose the 
conserved momenta Pi for the constant variables instead of the usual actions over a 
period § pdq that are more apt for bounded motions. Notice however that the Pi are 
different from the momenta appearing in perturbative calculations, even if both sets 
may coincide asymptotically or in some set of R". Coming back to our problem, the 
function W would be given explicitly as a complete integral of the following Hamilton- 
Jacobi equation: 

Due to the relations b. 1 and b.2 above, it is permitted to write W as the path-independent 
integral: 

n „q 

W{q,P)^Y. / dqiUq,C) (5) 
i=l ■''10 
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where go is a constant configuration space point, and the Ci (that remain fixed during 
the integration) are functions of the Pi whose determination is necessary to solve ex- 
pHcitly the problem. We are not concerned here with the search of specific solutions, 
but with the fact that integrability ensures their global existence. In fact, the equations 
(H) can be written in the form 
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Pa{q,P) ^ (kaiq.C), Qa{q,P)^^l dq, '^'''^^ ^ ' , a = l,...,n (6) 

i=l-'lo 

The first equation is simply the definition of the isolating integrals (b. 1). As a bonus, 
time can be given as a function of phase space in two alternative ways: either in terms 
of the old variables, or equally in terms of the new ones. Consider that a particle 
initially at (q, p) arrives at the position q(i) = x in the instant t(x) = t, then: 



t(x)^^(X.-Q.) = ^5: fd^.^^gLZ), a=l,...,n (7) 



where Xa = dW{x,P)/dPa (obviously, Xa = Qait{x.)) by construction). Note that 
in (^ there is no summation over the index a. In fact, integrability can be envisioned as 
the simultaneous existence of n independent flows each of them contained in a different 
phase space plane. The requirement of integrability was noticed by Einstein who 
analyzed its implications for the old Bohr-Sommerfeld quantization conditions, that he 
reformulated accordingly giving a new condition, that was criticized by Epstein 
Integrability [|[| allows n different expressions to define the unique time of arrival. 
Only a pair {Qa, Pa) appears in each of them, and they all are equivalent. This holds 
even when there is no separable solution to the original Hamilton- Jacobi equation (Q) 
due to the presence of the potential V{q) in the Hamiltonian. Only for some well 
known cases [^] the problem is separable in the original variables. Independently of 
this, notice that as (Q(t), P) defines a straight line in the phase space, it is simple to 
lay one of the axes (the n*^ say) along it. This amounts to define Ho{P) = P„ which 
gives Pn{t) = E and Qn{t) = t + Qn, while the other variables remain constant 
Qj{t) — Qj, Pj{t) — Pj, j — 1, ■ ■ ■ ,n — 1. With this choice, one can write: 

\ r A dpziqi---qn,Pi--.Pn-i,E) 

- 2^ / dq^ (8) 

1=1 -'1 

with the Pi's given in (^. This is the standard equation of time that appears in the lit- 
erature. The rest of the relations would give the time independent geometric properties 
of the trajectories. Note that for central potentials only depends on E, so that (^ 
reduces to t{rx) — j^^ dr{dpr/dE). 

We have focused the discussion of this section on the dual definition of the time of 
arrival, that can be given in terms of the original phase space variables, or of the free 
translation variables. This duality is a foundation stone for the quantum method pre- 
sented in this paper We will obtain the time of arrival operator of interacting particles 
t{x) by applying a quantum version of the canonical transformation W{q, P) to the 
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well known operator for the time of arrival of free particles io{x). The properties of 
the latter have been extensively analyzed in the literature. For completeness, and to fix 
the notation, we present a summary of them in the next section. 



2 Time of arrival of free quantum particles 

In one space dimension Eq. (Q) gives the time of arrival at a; of a free particle initially 
at {q,p) as a function of the phase space variables that depends on x parametrically: 
tQ{q,p;x) = m{x — q)/p. In spite of its simplicity, this expression presents serious 
quantization difficulties |Q, |[ |, [lO[ |ll|] whose solution we oudine here [|l^ |l2[ [l3[ |l^ 



|15p . First of all, it requires a decision about operator ordering, the simplest one being 
symmetrization: 

^o(g,p;a;)=m(^-i{g,4}+) = -e-'^^/^ qMe^^^ (9) 
p 2 p \ P \l P 

Notice the proliferation of carets above. It is a reminder that we now deal with operators 
acting on the Hilbert space of the free particle states. From now on, we will drop the 
operator carets, simplifying the notation as much as possible, wherever this will not 
produce confusion between operators and c-number variables. The eigenstates \txsQ) 
of this operator to{x) in the momentum representation can be given as (Ti ~ 1) 

{p\txsO) = 6i(sp)wM exp(i;^i) {p\x) (10) 
V m 2m 

where t is the time eigenvalue, x the arrival position, and where we use s ~ r for 
right-movers (p > 0), and s — I for left-movers [p < 0.) The label stands for free 
particle case. Finally, the argument sp of the step function that appears on the rhs is 
+p for s ~ r, and — p for s = /, so that 

/■oc 

Oirp)^ / dp\p){p\ and e{lp) ^ / dp\p){p\ (11) 

Jo J-oo 

The degeneracy of the energy with respect to the sign of the moment is explicitly shown 
by means of a label s — r,l in the energy representation, where 

{Es'0\txsO) = Ss's ( — )i/4e'^* {EsO\x) (12) 
m 

Summarizing, there is a representation for the time of arrival at x spanned by the 
eigenstates 

ItesO) = (^)i/4e.//o*n,o|^) (13) 

TO 

where IIso projects on the subspace of right-movers (s = r), or of left-movers (s = I), 
i.e. 

/•OO 

n,o = / dE\EsQ){EsO\=0{sp) (14) 
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These time eigenstates are not orthogonal. This gave rise in the past to serious 
doubts about their physical meaning. The origin of the problem can be traced back 
to the fact that is not self-adjoint, that is, that {(p\tQ{x)ip) ^ {tQ{x)ip\'ilj). This was 
pointed out by Pauli |Q] long time ago and is due to the lower bound on the energy spec- 
trum that prevents the applicability of the Stone theorem JT^. The problem emerges 
as soon as one attempts integration by parts in the energy representation. 

Not being self-adjoint or orthogonal, this operator poses an interpretation problem 
that can be solved by considering it in terms the Positive Operator Valued Measures 
(POVM). This is a class of operators less restrictive than the traditional projector valued 
measures. The POV measures only requires the hermiticity of to (a;) (i.e. to{x) = 
{to{x))*^) to assure the positivity of the measure. Now, instead of a Projector Valued 
spectral decomposition of the identity operator, one has the POV measure 

r2 

PoiU{x);ti,t2) = y2 dt\txsO) {txsO\ 

Ji m m 

whose notation indicates the arrival interval and that the dependence on the arrival 
position comes through the projector ll{x) ~ \x) {x\ on the position eigenstate. For 
the above measure Po{l, 2)^ ^ ^o(l, 2) because \txsO) {txsO\ is not a projector, as the 
states are not orthogonal. However, the limit as t ^ cxd of Po{—t, +t) is the identity as 
can be checked explicitly. The time operator obtained is well suited for interpretation. 
This solution was introduced in [Q, and extensively analyzed in refs. [|l7[ |l8]]. It has 
been recently reviewed in [ |l9| ] and criticized in [ ^o| ] . 

In this formulation the time of arrival is given by the first moment of the measure 

+ 00 



/-l-oo 
dtt \txsO) {txsO\ 
'OO 



+OC Q TT o rr 

dtt{fi^f/^ e'^o* Va{x) e-»^°* (±^££)i/4 (16) 

-oo "1 TO 

where Vq{x) = IIso n(a;) IIso, which is not a projector. We now have the tools 
necessary for the physical interpretation of the formalism: Given an arbitrary state ^ 
at < = 0, its time of arrival at a position x has to be, according to ([1^, 

(^lio(x)IV^) = — / dU|(tesO|^)p, (17) 

with the standard interpretation of | {txsQ\%l}) p like the (as yet unnormalized) prob- 
ability density that the state j-)/)) arrives at x in the time t. The probability of arriving at 
X at any time is then P{){x) = j dt J^s I {txsO\ip) p, giving a normalized probability 
density in times of arrival 

Po(i,x) = --^ 5]|(txsO|V)P (18) 
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normalization that has been used in (p7|). 

The above equations ( |l7| , p^ can be given forms that are very useful for compu- 
tation and that throw some light on the physical meaning of the different quantities 
involved. By using explicitly (O), one gets 



Poix) = y^{fdEi—)'^^{x\EsO){EsOm* 
^ — ^ / m 

s 

= 27ry / dE{ — )^^^\(x\EsO){EsO\ib)\^ (19) 
^-^ J m 



The use of a similar procedure in (17 ) leads to {■>p\to{x)\'ip) 

= -TTTtE / dE{—y/mx\EsO)(EsOm* ^ {{x\EsO) (EsOm (20) 

s 

= irr-^ E / dEi — nx\EsO){EsOM' ±.{e,rg{x\EsO) + arg(i?sO|V)} 
Fq[x) ^ J m oE 

This expression is easy to understand. In fact it involves two ingredients: the plane 
wave amplitude (a;|£'sO) — exp(ispa;), along with the bracket (EsO\ip) = 

}ip{sp) where ip is the Fourier transform of the initial state in momentum space, 



and p — \/2mE. This gives for the arrival amplitude 



1 I 

{txsO\ij)^— dp J-e'^-^'+'P='^{sp) (21) 
27r ./o V m 



This is a free case so that the probability of ever arriving to x has to be one. In fact 

POO r-\-oc 

Poix) = E / ^ / '^P \^^P^\' = 1 (22) 

where we used that in our notation rp ~ +p and Ip ~ —p. We also have: 

A{arg(.|^;.0) + arg(i^.O|V.)} ^^{sx- ^^^iMMl) (23) 
oE p op 

There are initial wave packets centered around the values (go,Po) for which ip{p) = 
\i^{p)\ exp(—zp(7o) with the amplitude peaked around po- Then, {d/dp) a.rgtp(sp) — 
sqo, and the time of arrival at x reduces to 



/"OO 

MxM = V / dp\i,isp)\ 



,m{x-qo)_ r^dp\m\'""^"~'"^ 



SP J-oc P 

(24) 
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which is the time of arrival of the classical free particle averaged over its initial state. 

The generalization to the case of three space dimensions [ [l3| | is not straightforward. 
The reason is that to begin with, there are three equivalent equations for the time of 
arrival. To be compatible, they have to satisfy the constraints 

£ = (q-x)Ap = (25) 

where we drop the distinctions made in (^ between upper and lower case letters, as 
they are the same objects for free particles. Classically, the constraints correspond to 
the fact that x has to be a point of the particle's trajectory, therefore the angular mo- 
mentum can be written as L = x A p. In other words, the angular momentum with 



respect to x, that is £, has to vanish. We now show that the constraints (25) are first 
class. First of all, they are closed as their components La = tahc {q — x)bPc satisfy 
the algebra of 3-D rotations, namely {Ca, Cb} = Cabc ^c- Then, the total Hamiltonian 

2 

is Ht — + X ■ C, where A is a vector multiplier, so that {Ca, Hx} = ^abc Xb Cc- 
Therefore, the constraints form a first class system that depends parametrically on x, 
one for each arrival position x. Not all the x's can be reached from an arbitrary set 
(q, p) of phase space variables. Only those x that satisfy the constraints are positions 
where the particles with these dynamical variables can eventually be detected. A de- 
tector placed somewhere else will miss them. 

The above translates into quantum mechanics as it is: Not all the states in the 
Hilbert space of free particle states H with Hamiltonian Hq can be detected at a spe- 
cific position X. Only the subspace Hx of the states that satisfy the constraints ( ^^ 
(where now q and p are operators) qualify as the Hilbert space of detected states (at x). 
This subspace is spanned by the states x) e 7i of the form x) = ip{Ho) \ x). 
Here, V'(^o) is an arbitrary function of Hq, that may also depend on other parame- 
ters, I x) is the eigenstate q|x) = x|x) of the arrival position. In particular, the 
detected subspace Hx is obtained from Ho by a translation of amount x, as required 
by covariance. 

The value of t comes from the equation of motion in the subspace orthogonal to 
the constraints, namely p • x = p • tQ{x)+cij, that in spherical coordinates where 
|p) = \p, 9p, (j)p) and q ~ i-^ can be written as x = ^ to{x) + q. This can be readily 
inverted to give 

to(x) = -- e-^P'^p-i/^ g^px (26) 

P 

Notice the characteristic powers of p to the right and to the left of q. This operator 
ordering makes of to a maximally symmetric operator with respect to the measure d^p, 
making integration by parts a straightforward task. In d space dimensions we would 
have to oc (~*^)-P" withn = (d — 2)/2 [[l3|]. The eigenfunctions of to are given 
in the momentum representation by: 



where t S R is the time eigenvalue, and Ep ~ p^ /2m. One can define a time of arrival 
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representation given by 

\t; X, 0) = -yL= (^^_) V^e'^^o * I ^) (28) 
V 47rm 2mHQ 

These eigenstates are not orthogonal. They correspond to a POV measure defined by 
the spectral decomposition 

/ + 00 
dt\t;x,0)(t;x,0\ (29) 
-oo 

It can be immediately seen that for any state e Tix, and for arbitrary momentum 

P 

(p|lx|V;x) = (p|V;x) VpeR^ (30) 

Therefore, the operator Ix is a decomposition of the identity within the subspace of de- 
tected states Tix- The fact that Ix < 1, so that the decomposition is uncompleted, is the 
quantum version of the classical case where only a part of the incoming particles will 
(reach and) be detected at x. From here it is clear that our formalism is finer than that 
provided by the so-called screen operators [^, that would describe the arrival at a two 
dimensional plane put across the particle trajectories. In fact, these screen operators 
would correspond to a coarse graining of the present formalism, whose interpretation 
is analyzed in some detail in [|l3ll . 

The time of arrival can be given through the first momentum of the POV measure 

io(x) = J dtt \t; X, 0)(t; x, | 

= f dtt {-^)^^^e'""*\^){^\{—^y/^e-'"''' (31) 

whose similarity with the 1-D case (|l6|) is evident, and can be used as a guide to get the 
average time of arrival an other quantities of interest that were worked on in one space 
dimension. 



3 The arrival of interacting particles 

In this section we want to determine the effect on the times of arrival of a position 
dependent interaction between the particle and the medium, that we describe by a po- 
tential energy V{q). For instance, we want to consider the case of a barrier placed 
between the detector and the initial state. We would put a detector at x (at the other 
side of the barrier), and prepare the initial state |^) of the particle at t = (at this 
side of the barrier). We would then record with a clock the time t when the detector 
clicks. Repeating this procedure with identically prepared initial states, we would get 
the probability distribution P{t, x) in times of arrival at x. This is the same procedure 
used for the free particle case, the differences coming from the presence of the potential 
energy V{q). 
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To find the quantum time of arrival we will use what we know from the classical 
case: There is a canonical transformation from the free (iJg = ^) translation variables 

2 

{Q,P) to the actual variables {q,p) of the interacting situation where H = ^ + 
V{q). Time can be given equally in any of these two versions and we did already 
quantize the free version to in the previous section. Now, in successive steps, we do 
the following [pi]]: We first construct the quantum canonical transformation U that 
connects the free-particle states to the eigenstates of the complete Hamiltonian. This is 
the quantum version of the (inverse of the) Jacobi-Lie canonical transformation (^|]). 
We will see later on that U is given by the Moller wave operator We will then apply U 
to to to define the time of arrival t in the presence of the interaction potential V{q) in 
terms of <o- We will work out the details of this transformation t = UtoW. Finally, 
we will also address some questions of interpretation of the resulting formalism. 

Dirac introduced canonical transformations in quantum mechanics in a number of 
different places by means of unitary transformations U (UU^ ~ U^U = 1 ). To 
fix the notation, we assume in what follows that the operators q and p are given in the 
coordinate representation of the Hilbert space L^{x) by q = x and p = —ih-^. If the 
operators q and p are the result of an arbitrary canonical transformation applied to q 
and p, then there is a unitary transformation U such that 

q^UUU, p=U''pU^ [q^ p] ^ [q^ p] = -ih (32) 

One can also define implicitly the quantum canonical transformations as is done in 
classical mechanics, a possibility that has been thoroughly analyzed and developed. 
The main results of the method are collected in [^^, which also includes references to 
other relevant literature. The definition of U is given implicitly by the two conditions 

F{q, p) = Fo(q, p), and G(g, p) = Go(q, p) (33) 

where F, G, and Go are functions of the operators shown explicitly as their argu- 
ments. They can not be chosen arbitrarily, the necessary and sufficient condition for 
the canonicity of the transformation being [F, G] = [Fq, Gq]. The dependence of (^3|) 
on U can be explicitly given by using (|2|) in it: 

F{q,p)U ^ Fo{q,p), ^MU^ G{q,p)U = Go{q,p) (34) 

that comes from the straight application of ( ^2[ ) to the first members. In addition, U is 
unitary so that the spectra of the original and transformed operators have to coincide. 
We now assume that F and F^ are self-adjoint operators whose eigenvalue problems 
are solved by the states |/ s) and | / s 0) (both corresponding to the same eigenvalue), 
that form orthogonal and complete bases of the Hilbert space satisfying 

F\fs)=\f\fs), Fo|/sO) = A/l/sO) (35) 

We are accepting here the presence of degeneracy indicated by the discrete index s, 
something that we will need later Assuming now a continuous spectrum (the case we 
will be interested in), the operator U that satisfies the first row of (|5]) is given by 

C/ = V/ d\f\fs){fsO\ (36) 
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It is straightforward to verify that it is unitary. We can now give the definition of G in 
terms of Go using U, that is, G = U GqU\ which in full detail reads 

G{q,p)^y^[ dXfdXf,\fs){fsO\Go{q,p)\fs'0){f's'\ (37) 

This is the main result of our procedure. The fact that we can define an operator G, 
canonically conjugate to F, if we know Go and U . 

We will now apply this to the case where Fq is the free Hamiltonian i/o, F the 
complete Hamiltonian H and Go the time of arrival of the free particle Eqs.(H), or 
([l|). Then, we have Hq^U^ HU and Ho (a;) = f/t n(x) U. Associated to the free 
particle there was the positive operator valued measure Po of Eq.(p3|). Accordingly, the 
POV measure P of the interacting case will be given by (c/(^4|)) 

P(n(.T);ti,t2) = UP^{Ilo{x)-hM)U^- (38) 

then, the time of arrival operator in the presence of interactions (the G of (^7|)) is given 
by 

t{H, n{x)) = C/io(-ffo, no(a;))f/t. (39) 

We noticed above that the spectra of the original and transformed operators had to 
coincide. Now, a'{Ho) — R+ so that not all the Hamiltonians can be obtained from 
Hq by this procedure. In general, some fixing will be required to make the spectra 
coincide. Here we will only consider well behaved potentials (V{q) > Vq G R), 
vanishing appropriately at the spatial infinity. This ensures the required coincidence of 
the spectra, but introduces two solutions for U due to the existence of two independent 
sets of eigenstates of H: 



poo 

Ui±) - / dE\Es{±)){EsO\ = r!(±) 
„ Jo 



(40) 



these are the MoUer operators connecting the free particle states to the bound and scat- 
tering states. In presence the presence of bound states these operators would only be 
isometric, because the correspondence between eigenstates of H and free states could 
not be one to one. In our case V{q) > 0, there is one free state for each scattering state 
and conversely. Thence, the Moller operators are unitary. In this case, the intertwining 
relations Hfti^j--^ = il(--j--)iJo can be put in the more desirable form H — rij^-ji/o^J^-j- 
We will also follow the standard sign conventions, choosing fi(+) in ( p^ that, when 
E = limg^o+ {F + gives signal propagation forward in time. The results that 
would be obtained with f^(-) would correspond to the time reversal of this situation. If 

r is the time reversal operator, then (n(x); ti, = -P(+)(n(a;); — ^2, ^ ^i) 
For notation simplicity, we will omit these labels (±) wherever possible. 



The parameter x that appears in (^7|) and (39 ) is the actual detection position in the 
interacting case, the place whose time of arrival at we want to know. Therefore, the 
arguments of t in ( ^9| ) have to be n(x) — \x)(x\ and H. Hence, the argument of to 
will be an object Ho (a;) ~ ri^n(x)il which collects all the states of the free particle 
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that add up to produce the actual position eigenstate \x) by the canonical transforma- 
tion. Much of the difference between the classic and quantum cases is hidden here. 
In particular, the quantum capability to undergo classically forbidden jumps in phase 
space has much to do with the fact that n(a;) and llo{x) can not be position projection 
operators simultaneously. 

We have now at hand all the tools necessary to answer the questions about the time 
of arrival of interacting particles. Given a particle that was initially (at t = 0) pre- 
pared in the state we can compute the predictions for the average time of arrival 
(V'|i(a;)|V'), the probability distribution in times of arrival P{t, x) and the probability 
of ever arriving at x, P{x). Instead of writing more equations, we refer the reader to 
Eqs. ( [I6HT7I ), ( |l5| , p^ and (|l9|). By simply erasing the label from them, one gets 
the correct expressions for the interacting case, with the caveat that -to be of practi- 
cal use- they require the knowledge of the scattering states and Moller operator It is 
worth to recall here that the expression ( pO[ ) for the average time remains valid after 
dropping the O's. So, {%lj\t{x)\%l}) is still the sum of two independent pieces, one con- 
taining [d/dE] aig{E s | V') that only depends of the initial state, the other that contains 
(d/dE) aTg{x \E s ) and only depends of the position of arrival. 

We now consider the case where there is a finite potential energy starting at the 
origin iV{q) = 0, Vq < 0), which is so smooth that the quasi-classical approximation 
is valid. Then for E > V{x) the exponentially small reflection amplitude can be 
neglected, giving the scattering states 

+ ,41, 



with p{q) = ^2m{E — V{q)), that are normalized to an incoming right-moving par- 
ticle by time unit. We now consider the physically interesting case where the initial 
wave packet is normalized to 1, (i.e. that / dp\4>{p)\'^ — 1 with ijj{p) — also, 
that it is localized around a position qq well to the left of the origin, and that it has a 
mean momentum po 3> V{x). Then, to this order the probability of ever arriving at x 
(c.f. (|T^) gives 

P{x) « e{-x) + e{x) P>{x), where Py{x) = / dp-^ \^Ap)? (42) 

Jo P{x) 

so that is the (unnormalized) probability of arrival at the point x with 

momentum p{x), as corresponds to the quasiclassical case. Notice that to the left of 
the origin the result is the same that in the free case. This comes about because the 
approximation neglects reflection, thus missing at g < any information about a the 
existence of a finite V al q>Q. For the time probability distribution one gets 



9ix) 







2'kP>{x) 







V p{x) 



(43) 
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which, not surprisingly, is the same as that of free particles for a; < 0. Finally, 

f°° ~ m 
i^PlU^) « 0{~x) dp \^{p)\^ -{x ~ qo} 
Jo P 

0{x) p I ~ i^^a f "^90 , _ dq 



dp-^ {-^ + m / ^} (44) 

P>(x) Jo P(x) P Jo PW 

Therefore, we recover the time of arrival of the free particles for negative x. On the 
other hand, for a; > we get the classical time of arrival at a; for initial conditions 
{qo,p), {m/p{q)) dq, weighted by the probability of these conditions. 



4 Advanced or delayed arrival? 

What is the effect of putting a quantum barrier in the path of the arriving particle? Hart- 
man studied this question a long time ago, reaching the conclusion that tunneled 
particles should appear instantaneously on the other side of the barrier Our formalism 
supports this result, but only for thin enough barriers. 

The time of arrival at a point x in the presence of a barrier will be given through a 
probability amplitude 

(^51^-)= f dE{ — )^/^e-''''{x\Es{+)){Es{+)\^lj) (45) 
J m 

In the case where x is at the right of the barrier, the amplitude can be approximately 
given by Q 

(tes|V>) « 4^ / dE{^y/^e~^(^'~P^^Tip)^jip) (46) 

where T{p) is the transmission amplitude for momentum p. Now, the total probability 
of eventually arriving at x in any time t is P{x) ^ J dp \T{p) ■tp{p)\'^ that is inde- 
pendent of X in cases like this, where x is beyond the range of the potential. After a 
straightforward calculation we get for the average time of arrival at the other side of 
the barrier the corresponding version of ( ^o| ) 

(^|t.|V')«7jK r ^iP)\"-{^ - go + ^"''gf ^^^^ } (47) 

P{x} Jo p dp 

It is the value of the Wigner time [ p6| ] averaged over the transmitted state. 

Consider a simple square barrier of height V and width a. The transmission coeffi- 
cient is in this case: 

T{p) = e-'P'' (i_^ (P^+P'^^ temp' a] secp'a (48) 
V 2pp' J 



where p' = ^/p^ — Py, that is imaginary for p below pv- Notice the contribution —pa 
to arg(r(p)). This will subtract a term a to the path length x — qo that appeai-s in 
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The barrier has effective zero width or, in other words, it is traversed instantaneously. 
This is the Hartman effect for barriers. To be precise, the effect is not complete, it is 
compensated by the other dependences in p'a present in the phase of T{p). In fact, 
it disappears for low barriers (pv/p) 0, where all the a dependences of the phase 
cancel out, as was to be expected because the barrier effectively vanishes in this limit. 
In the opposite case of high barriers {p/pv) the effect saturates and there is a 
decrease — ^ in the time of arrival of transmitted plane waves, that emerge almost 
instantaneously at the other side of the barrier 

The averaging of the Wigner time over the transmitted state, present in (^7[) as a 
consequence of the formalism, has dramatic effects, because it effectively forbids the 
transmission of the wave components with low momenta. In fact, it produces the expo- 
nential suppression (by a factor exp(— 2|p'|a)) of the tunneled components. Therefore, 
only the components with momentum above pv have a chance of surmounting thick 
barriers, being finally transmitted. But these components are delayed by the barrier 
(for them ( "^^'^^J^^-*-* ) > 0), whose overall effect transmutes from advancement into 
retardation [21 1 at a definite predictable thickness that depends on the barrier height 
and also on the properties of the incoming state. 
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